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We consider an array of coupled cavities with staggered inter-cavity couplings, where each cavity mode 
interacts with an atom. In contrast to large-size arrays with uniform-hopping rates where the atomic dynamics 
is known to be frozen in the strong-hopping regime, we show that resonant atom-field dynamics with significant 
energy exchange can occur in the case of staggered hopping rates even in the thermodynamic limit. This effect 
arises from the joint emergence of an energy gap in the free photonic dispersion relation and a discrete frequency 
at the gap's center. The latter corresponds to a bound normal mode stemming solely from the finiteness of 
the array length. Depending on which cavity is excited, either the atomic dynamics is frozen or a Jaynes- 
Cummings-like energy exchange is triggered between the bound photonic mode and its atomic analogue. As 
these phenomena are effective with any number of cavities, they are prone to be experimentally observed even 
in small-size arrays. 
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I. INTRODUCTION 

The tremendous interest spread over the last few years in 
the dynamics of coupled-cavity arrays (CCAs) |Q]] has made 
this emerging field a topical one within the general framework 
of quantum coherent phenomena and beyond. Many are the 
reasons behind such a substantial and widespread attention. 
Along with cold atoms in optical lattices CCAs stand out 
as an attractive controllable test bed for many-body phenom- 
ena such as quantum phase transitions ^ . An appealing fea- 
ture of CCAs is local addressing, namely the available high 
control of each single site in terms of performable measure- 
ments, quantum-state engineering and dynamical parameters 
tuning. From the perspective of quantum optics, CCAs open 
the door to the exploration of regimes where an almost ubiq- 
uitous feature of so-far-investigated atom-photon dynamics, 
i.e. the effectiveness of descriptions in terms of single-atom 
dynamics, does not hold any more. Importantly, the variety 
of possible experimental setups prone to implement CCAs |4|] 
and, mainly, the widespread expectation that the technology 
required for their actual fabrication is by now at hand (at least 
for small clusters) are providing formidable motivations to ex- 
plore the rich physics of CCAs. In particular, one of the lines 
along which current investigations are proceeding is the study 
of excitation propagation in various forms such as transport of 
photons and excitons SB!,polaritons @,[T1 and solitons 
ifliri . In this framework, transport in CCAs with defects or im- 
purity atoms is also receiving considerable attention lfl2l - [l5ll . 

In a typical arrangement of CCAs, the field mode of each 
cavity is coupled to a two-level atom via a Jaynes-Cummings- 
type interaction. Energy can be thus stored in the form of 
both photons and excitons (atomic excitations) that are able 
to transform into each other as well as propagate along the 
array. While photons can move by direct hopping between 
neairest-neighbour cavities, excitons can travel only provided 
that they are transformed into photonic excitations and even- 
tually converted back. The dynamics ruling excitation trans- 
port in CCAs is therefore non-trivial and physically attractive. 
While first works along this line addressed small clusters jij.^] 



only more recently arbitrary-size arrays have been tackled @|. 
In all such scenarios, for uniform values of atom-photon inter- 
action strength, cavity-mode frequency and atomic detuning 
and whenever the decoupling of the field's hopping Hamil- 
tonian in terms of normal modes is known the Hamiltonian 
describing the full many-body dynamics enjoys an attractive 
feature. It can indeed be rearranged jj as the sum of de- 
coupled Jaynes-Cummings (JC) models uM, each coupling a 
field (bosonic) normal mode to its excitonic analogue. So long 
as the array size is small (assuming non-degenerate photonic 
normal modes) resonant excitation of only one of such effec- 
tive JC systems is possible by judicious tuning of the atom 
frequency 0]. Physically, this circumstance is quite remark- 
able since it entails that significant amounts of energy can be 
exchanged between atoms and photons even in the strong- 
hopping regime, i.e. when the atom-photon interaction rate 
is much lower than the photon hopping rate. Such a pic- 
ture, however, can drastically change with large-size CCAs. 
By taking uniform-hopping rates, for instance, the passage 
to a large number of cavities brings about that the free field 
normal-mode spectrum tends to a continuous band. Clearly, 
this rules out the possibility of any resonance and thus first- 
order atom-field energy exchange. As a result, in the regime 
of strong hopping the atomic dynamics turns out to be frozen 
JH]. One may wonder whether the above picture still holds 
for large-size arrays when a non-uniform pattern of hopping 
strengths is considered. Motivated by some findings in the 
framework of spin chains [17] Makin et al. ^ investigated 
a parabolic distribution of hopping rates and found no signif- 
icant changes in the strong-hopping regime compared to the 
uniform-hopping behavior. 

In the present work, we show that resonant atom-photon in- 
teraction can take place even under strong-hopping conditions 
in an arbitrary-size array with a staggered pattern of hopping 
rates. In particular, we demonstrate that when certain array 
sites are initially excited and each atom has negligible detun- 
ing from the single-cavity field a JC-like dynamics involving 
a bound photonic normal mode interacting with its excitonic 
analogue is triggered. Among other effects, it entails a signifi- 
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cant exchange of energy between field and atoms. This behav- 
ior basically stems from two features. First, the gapped nature 
of the normal-frequency spectrum of the hopping Hamilto- 
nian. Second, the occurrence under open boundary conditions 
(BCs) of a discrete frequency corresponding to a bound nor- 
mal mode at the center of the gap. Such circumstances make 
resonant atom-photon dynamics possible even in the thermo- 
dynamic limit. 

The present paper is organized as follows. In Section[ll] we 
introduce the set-up which we focus on, an array with stag- 
gered hopping rates, and give the Hamiltonian that describes 
its dynamics. We then focus on the special case where the free 
field Hamiltonian reduces to that associated with a uniform- 
hopping array, which allows us to briefly review the dynam- 
ics of such arrays. In Section [HI] we focus on the free field 
Hamiltonian in the general case, present exact solutions for 
its normal modes and associated frequencies and discuss its 
main properties. We then illustrate the regime which we fo- 
cus on and give the corresponding effective representation of 
the full Hamiltonian. In Section HV1 we show how an initial 
atomic excitation propagates along the array and discuss the 
salient features of the excitation-transport dynamics. Finally, 
in Section |Vl we give some comments and draw our conclu- 
sions. 



H. HAMILTONIAN AND REVIEW OF THE 
UNIFORM-HOPPING SETTING 

We consider a finite-length array of N low-loss cavities, 
where nearest-neighbour cavities are coupled together so as 
to allow for photon hopping. Each cavity sustains a single 
field mode of frequency a)/, which is coupled at a rate / to a 
two-level atom of Bohr frequency a> a (here and throughout we 
use units such that Ti — 1). The Hamiltonian of the full system 
reads 



H - Hf + H a + H, 



(1) 



where 



N-l 



H f = w/^ °% - «Yj [1 " + h ' c -) ' (2) 

(3) 

Hi = J^(h x &l + il& x ) . (4) 



N 

H„ = u) a ^b\b x , 

x=i 

N 



i.e. [b x , b x >] — [b x , b x ,] = for any x + x' since operators acting 
on different Hilbert spaces commute. On the other hand, it is 
easily checked that Vx hi = [b\] 2 =0 and b x b\+b\b x = 1. 

The distinctive feature of Hamiltonian ([1]) is the staggered 
pattern of hopping strengths. The odd (even) hopping rates, 
i.e. those associated with the nearest-neighbour cavities 1-2, 
3-4,... (2-3, 4-5,...), take the value K\ = (1 +rf)K [k 2 = (1 -tj)k)], 
where k is a hopping rate. Hence, two different hopping 
strengths are periodically interspersed along the array. As 
by setting r\ = an array with uniform hopping rate k is re- 
trieved [ 8] the dimensionless staggering parameter rj measures 
the distortion of such uniform-hoppings CCA. For \t]\ — 1, the 
array reduces to a collection of independent two-cavity blocks 
(but one comprising a single cavity in the case of odd N). A 
sketch of the whole set-up is given in Fig. 1(a). The distorted 
tight-binding model [Qjl specified by Hf has recently been 
harnessed (in its fermionic version) in Ref. lfl9ll for the sake 
of quantum state transfer. As we thoroughly discuss in the 
next section, regardless of the value taken by T], Hf in Eq. (f2]i 
can be exactly arranged in a diagonal form in terms of normal- 
mode field operators. In the remainder of this Section, though, 
we focus on the uniform-hopping case 77 = 0. 

For 77 = 0, the free field Hamiltonian can be expressed in 
terms of normal modes as Il8lll7ll 



where 



k = 



2nm 
N + 1 



(Ok a k a k 



(m = l,...,N) 



a>k — o)j + 2k cos — , 



a k - 



sin — x] a. 



(5) 

(6) 

(7) 

(8) 



In Eqs. (O and ©, a k and or, are, respectively, bosonic an- 
nihilation and creation field operators associated with the kth 
photonic normal mode, which satisfy standard commutation 
rules. 

A feature of Hamiltonian (Q]) is that atom-photon interac- 
tion strengths, cavity-mode and atomic frequencies are uni- 
form throughout the cavity array. This allows to arrange it 
in an elegant and useful form Jj, IM] in terms of decoupled 
effective JC models according to 



H 



O) k &l& k + (o a 0l$ k + j{&l$ k +0la k )] , (9) 



In Eqs. ©-(HI, a x (a x ) is a bosonic field operator that an- 
nihilates (creates) a photon at cavity x, whereas b x (b x ) is 
an atomic operator that annihilates (creates) an exciton on 

the xth atom according to b x = |/SiJ = \g) x {e\, where \g) x 
(\e) x ) is the ground (excited) state of the xth atom. Hf and 
H a are the free Hamiltonians of the field and atoms, respec- 
tively, while Hi describes the atom-field coupling. Notice 
that atomic operators associated with different sites commute, 



where the expansion of the atomic normal-mode operators 
/Vs in terms of site operators b x 's is apart from their differ- 
ent commutation rules fully analogous to Eq. ([8]). Notice that, 
unlike & k 's, each f} k has the same associated normal-mode fre- 
quency co a . 

In each effective JC-like Hamiltonian appearing in the de- 
composition (0, a field normal mode couples to its atomic 
analogue through a JC-type interaction. As all of the atomic 
modes have the same frequency a> a when co a = for a given 
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k' only the corresponding pair of field and atomic normal 
modes is, in principle, resonantly excited. 

For instance, for N — 2 only two normal modes arise in 
the decomposition (0. Their associated annihilation opera- 
tors and energies are in such a simple case conveniently ex- 
pressed as 0] a ± — (a i ±ai)j V2 and cj ± = (Oj±k, respec- 
tively [see Fig. 1(b)]. Analogously, the atomic normal oper- 
ators read fi± = {b\ ±1)2)1 y2. When the atoms are tuned in 
resonance with one of these two modes by setting A = +k 
(where A = a>f — co a is the atoms' detuning from the single- 
cavity frequency) a JC-like dynamics occurs with a continu- 
ous exchange of energy between the involved photonic mode 
described by a ± and the corresponding exci tonic mode 17J . A 
typical behavior fl] that arises under such conditions is that 
when one atom, say atom 1, is initially prepared in the excited 
state (with atom 2 in the ground state and no photons popu- 
lating the cavities) the photonic normal mode is progressively 
excited while the excitation probability of atom 2 grows. Over 
the following stage, while the field returns the received energy 
atom 2 keeps increasing its excitation probability until at a cer- 
tain time the initial excitation has fully transferred to it (with 
the field and atom 1 unexcited). Afterward, the phenomenon 
is reversed through a further excitation-deexcitation cycle of 
the field mode until the initial state is fully retrieved. Trigger- 
ing such a dynamics, which is of first-order in J, through a 
judicious tuning of the atomic frequency is however possible 
only in small-size arrays. 

Indeed, as N grows [see Fig. 1(b)] more and more photonic 
levels gather in the range [cjj — 2k, cjf + 2k] as implied by 
Eq. (0. In the thermodynamic limit N — > oo, the photonic 
spectrum takes the form of a continuos band. For large arrays 
this in fact suppresses any possible resonant atom-photon cou- 
pling due to the lack of discrete frequencies. As a result, in the 
strong-hopping regime J «: k a straightforward use of the in- 
teraction picture shows that so long as second-order processes 
are negligible the atomic dynamics becomes fully frozen @|. 
Hence, in striking contrast to the above-described process in 
a two-cavity array, an exciton initially localized on a given 
atom is unable to move away. Also, energy exchange between 
atoms and field is suppressed. 



III. NORMAL MODES FOR STAGGERED HOPPING 
RATES AND EFFECTIVE HAMILTONIAN 



The scenario discussed above for a uniform-hopping array 
substantially changes when rj + 0. To illustrate this, in Fig. 1 (c) 
we have numerically computed the normal-mode frequencies 
of the free field Hamiltonian © for t] = 0.5 and N = 3, 5, 15, 5 1. 
Two major differences appear compared to the case i] — 0. 
First, the spectrum is gapped and thus in the limit of large N 
two continuous bands emerge, instead of a single one. Second, 
regardless of N a discrete normal frequency lies exactly at the 
center of the gap. These findings can be made analytically 
rigorous. 

To begin with, for the sake of notation compactness we de- 



fine the two quantities 



7] + 1 



2a: A /cos 2 -+77 2 sin 2 - . 
2 ' 2 



(10) 
(ID 



Note that £j is defined positive. 

It can be checked that for an arbitrary odd number of cav- 
ities N (later on we comment on this assumption) Hj can be 
exactly arranged in the diagonal form 



(12) 



k ji=± 



where k = 2nm/(N + 1 ) with m = 1 , 2, . . ., (N- 1 )/2 and 



= — 2 V t* 
77-I Vt™-! jr{ 



07 

CO k ± = (Of + E k 



ak± = 



N+l 



^ sin (kx) a2x ±^] sm (kx+ftk) a.2x 



X=\ .1=1 

The phase §u appearing in Eq. ( fTBT ) obeys 



Kjl-T]) 

£k 



(13) 
(14) 

(15) 



(16) 



The discrete spatial functions specifying the expansion of 
normal operators Eqs. ([TBI and ( fTBI l in terms of site opera- 
tors {a x } can be shown to fulfill orthonormality conditions, a 
proof that we carry out in Appendix iBl As such, [at, &k,±) 
form a set of bosonic annihilation operators fulfilling standard 
bosonic commutation rules. As is to be expected, the normal- 
frequency spectrum {ojf,a)k±} and corresponding set of nor- 
mal annihilation operators {&£, &k±] reduce to (|7|i and (0, re- 
spectively, in the special case rj = 0. This check is carried out 
in detail in Appendix C. 

At variance with the problem tackled in Ref. lfl9ll that al- 
lowed assumption of cyclic BCs, in deriving Eqs. (fT2l-([T4li 
we have used open BCs, which are the natural ones to im- 
pose in a finite-length CCAs Jit]. Similarly to the uniform- 
hopping case lfl7ll . the solutions in the continuous spectrum 
[see Eq. $15[ ] are straightforwardly obtained by superposing 
those fulfilling cyclic BCs in the case of A^+ 1 sites lfl9ll and 
then requiring that they obey hard-wall BCs. The discrete so- 
lution &( is readily found by direct demonstration once one 
imposes that its normal frequency be a>f (so that the hopping 
term in Hf vanishes). In Appendix iDl we show in detail that 
upon replacement of ( [T4b . ( [TBI and ( [TBI ) on the right-hand side 
of Eq. (T% the free field Hamiltonian in Eq. (f2]i is retrieved. 

The two aforementioned bands in the normal-frequency 
spectrum are analytically described by Eqs. ( fTTT i and dT4b [see 
Fig. 1(c)] with their associated annihilation operators given by 
Eq. ( fT5l ). Their energy gap Aw, which coincides with the one 
obtained with cyclic BCs [HI], fulfills 



Aw > 4k\t]\ , 



(17) 



(a) 



(b) 



(c) 
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FIG. 1 : (Color online) (a) Schematic sketch of an array of N coupled cavities with staggered hopping strengths K[ = ( 1 + rj)K and k% = (1 - 77)*-. 
The field mode of each cavity is coupled to an atom (red dot), (b) 77 = (uniform-hopping case): normal-mode-frequency spectrum of the free 
field Hamiltonian Hf for N = 2,5, 15,51. Here and in next plot (c) frequencies are expressed in units of k and measured from <±> = u>f. (c) 
77 = 0.5 : normal-mode-frequency spectrum of Hf using open BCs for N=3, 5, 15,51. The levels below (above) the /c-axis [red (blue) dots] fall 
within the range [a>f—lK, iOf—K\ ([u>f+K, W/ + 2/C]), while a single discrete frequency a> = a>f (green dot) lies at the center of the gap between 
such intervals (its associated abscissa has been arbitrarily assigned). 
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FIG. 2: (Color online) Amplitude of the bound mode in Eq. J 1 3 1 > 
against the cavity number for N = 51 and 77 = 0.2 (a), 77 = -0.05 (b), 
77 = -0.2 (c) and 77 = -0.5 (d). As the array length N only affects the 
normalization factor, when N is varied the shape of these curves is 
unaffected. Notice that the amplitude vanishes at even sites. 



where the identity occurs in the thermodynamic limit N— > 00. 

For 77 — > the gap width vanishes so as to retrieve the 
uniform-hopping case ^ analyzed in the previous section. 
As \tj\ grows the lower bound of the gap on the right-hand side 
of ( [P7| i linearly increases at a rate proportional to k. As for the 
normal mode in Eq. ( fT3] l, its associated frequency u> j lies out 
of the two bands. Indeed, such normal mode is bound in that, 
differently from modes (TT3T >. it does not extend over the entire 
array. This is evident from the exponential functional form 
in Eq. ( fT3l and Fig. 2, which shows that the mode amplitude 
decreases (increases) exponentially with the cavity number x 



for negative (positive) values of 77. Notice that while the mode 
amplitude strongly depends on 77 its frequency is fully inde- 
pendent of both 77 and k. Using Eqs. ([Tol l and ( fT"3l l, the charac- 
teristic length A over which the bound mode is spread is found 
as A — 1 1 / log[( 1 +?7)/(l — /7)]| . Thus the local-mode length be- 
comes infinite when 77 0, as expected, while in the range 
< [77I < 1 it progressively shrinks for growing I77I until it be- 
comes negligible when I77I 1. This behavior is highlighted 
in Fig. 2, where the mode's spatial profile for different 77's is 
reported. A further feature which is evident in Eq. ( fT3l l and 
Fig. 2 is that regardless of 77 the bound mode vanishes at even 
sites. 

We now turn our attention to the full Hamiltonian ([TJ in 
order to highlight the implications of Eqs. ([T2"li-([T5ll on the 
CCA's dynamics. We first define atomic normal operators in 
full analogy with Eqs. (fT3l and (fT5l l as 



h± -- 



2 

N+l 



2^ sin (kx) l>2x ± 2^ sin (fcc+^it) &2x- 



(18) 



(19) 



where again k — 2nm/(N+ 1) with m = 1,2, (N— 1)/2 with 
and r still given by Eqs. ( fTol l and ( [Tol l, respectively. It 
is worth pointing out that unlike {df.Qty} the set 0c,^k/j) do 
not obey commutation rules since, while atomic site operators 
associated with different cavities commute, b x b[ — 1 - b\b x 
for any x. In analogy with the uniform-hopping case in the 
previous Section, the full Hamiltonian in terms of field and 
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atomic operators {&(, &k±,pt,Pk±) takes the form 



IV. DYNAMICS OF EXCITATION TRANSPORT 



H = H f +^[w t± a^ff^+w flJ 6^^ + y(j6^a^+h.c.)| , (20) 



kfj=t 



where 



ti ( = a> f a]ae + a> a jP ( Pt + J {p\a c +a e p\) . (21) 



As we show in Appendix [El the full Hamiltonian normal- 
mode decomposition < f20b straightforwardly follows from de- 
composition ( fT2l and the fact that to a and J do not depend on 
the cavity site. 

We take as free Hamiltonian Hq = Hf + H a [see Eqs. (f2]i 
and (0)]. Due to [b x , oy] = for any x + x', the commuta- 
tor between each site atomic operator and the free Hamilto- 
nian is given by [b x , Hq] = u a b x , which immediately yields 
[/^, Ho] = tJafif (£ = (, {k,[i}). Hence, in the interaction pic- 
ture each normal atomic operator evolves with time according 
to 0J\t) = fee-™"' (£ = £,{k,n\). On the other hand, using 
Eq. ( U21 the normal field operators evolve in the same picture 
as af'it) — ate~ ,0J f' and a^{t) — &k M e~ U0 'v t . In the interaction 
picture, the interaction Hamiltonian [cf. Eq. d20"l )l thus reads 

Hf = J&l [} { e iAt + {&l K <***-**) + h.c. (22) 

k/i=± 

When A = 0, i.e. each atom is on resonance with its own cavity 
mode, Eq. ( 122b shows that while the contribution to Hi from 
the bound mode is constant in time, those arising from the 
normal modes in the two energy bands rotate at frequencies 
that are at least equal to half the energy-gap lower bound 2k\t]\ 
[see (ITTbl - Hence, provided that 



J/k «: \t}\ , 



(23) 



such contributions are rapidly rotating and thus do not affect 
the system's dynamics. When A ^0, this conclusion still holds 
provided that the detuning is at most of the same order of mag- 
nitude of J and obeys |A| <k k\t}\. Therefore, in the above 
regime the effective Hamiltonian in the Schrodinger picture 
reduces to 



H 



eff 



H e + ^ j {coki&l^a^ + 1 

k/j=t 



(24) 



Eq. d24b embodies a central finding of this work, namely 
the possibility that the complex many-body atom-photon in- 
teraction reduces to an effective JC-like coupling described by 
H( between a bound photon mode and its excitonic analogue 
with all of the remaining field and atoms' modes freely evolv- 
ing. According to Eq. d23l ), when \r\\ < 1 such as in the nu- 
merical examples in Fig. 2 the occurrence of this few-degree- 
of-freedom dynamics is fully compatible with the strong- 
hopping regime J <sz k. In the thermodynamic limit, i.e. in 
practice for large-size CCAs, this marks a major difference 
from the uniform-hopping array {U discussed in the previous 
section. 



Our next goal is to shed light on the features of the sys- 
tem's time evolution in the one-excitation subspace in line 
with Refs. 0Htl- We recall that, despite normal atomic oper- 
ators associated with different modes in general do not com- 
mute, in the one-excitation subspace they effectively do 12011 . 

We denote by |0) the system's state with zero excitations, 
either photonic or atomic, and define states |*P±) as 



\y ± )=A ± a'' ( \0) + B ± pl\0) , 



with 



2.1 



B + = 



V(A ± Q) 2 + 4J 2 
A±D 

V(A + Q) 2 + A J 2 



(25) 

(26) 
(27) 



Here, Q = VA 2 + 4 J 2 is the usual Rabi frequency associated 
with the JC-like Hamiltonian (|2TT i. \*¥ ± ) are eigenstates of Hf, 
and hence of H e ff, with eigenvalues (u> a +u)f)/2 ± £2/2. Using 
this and taking as paradigmatic initial state ^(O)} one such 
that a given atom xq is excited, i.e. ^(O)) = b' Xo |0), at a later 
time t the system has evolved according to 

mt)) = <0|yS f ^J0>(B +e -''" /2 '| > I' + >+B-e'' n/2 '|^->) 

+ XX < 01 ^ fi * |0) ^^^L |0) ' (28) 

k n=± 

up to an irrelevant phase factor. Eq. (|28| | was obtained by 
expanding ^(O)) in the basis of stationary states of H e ff 
{\ x ¥±) , &k,± |0) ,Pk,± |0)}. In the following analysis of the im- 
plications of Eq. d28l >. we shall make use of the completeness 
of the basis of single-photon states {or! |0) , a, |0}}, i.e. 



k H=± 



I Olkj! = 1 lph 



(29) 



where 1 i p h is the identity operator in the one-photon Hilbert 
space of the field. Eq. d29b straightforwardly follows from 
orthonormality identities dB7b - dB9b . 

Let us first consider the case that xq is even. As the bound 
mode ( fT3l does not overlap even sites (see previous section) 
we trivially have <0| ySf bL |0) = so that only terms propor- 
tional to P k |0) do contribute to I'Fff)). This immediately 
yields that (0| a x \¥(t)) = for any x, i.e. no field excitation 
is developed. As for <0| b x ^(f)}, namely the probability am- 
plitude to find the xth atom excited, use of Eq. ( [291 , which is 
clearly valid for operators yS^ ± 's as well, along with Eq. ( TTBT ) 
entails (0| b x |*P(f)) = £>a,.i - ^ n °th er words, when xq is even 
PF(f)) — ^(O)), i.e. the atomic excitation is frozen analogously 
to the uniform-hopping case H]. The freezing behavior how- 
ever may not occur when xq is odd. Indeed, through a rea- 
soning similar to the one carried out above it is immediate 
to prove that for odd Xq if x is even then both (0| a x |*P(f)) 
and (0| b x \¥(f)) vanish, namely the initial excitation can only 
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spread over odd cavities. On the other hand, when x is odd 
projection of Eq. ( l28l i onto a x |0) and h\ |0) respectively yield 



\a x m)) = -2iNr^^ sin "r <T 



n a n 

cos— /-z— sin— f 
2 Q 2 



1 



(30) 



,(31) 



where N is the square of the factor out of the sum in Eq. ( fT3l ). 
Eqs. ( f30b and QTb fully describe the time evolution of the ar- 
ray in the regime such that the system Hamiltonian is well 
approximated by H e ff as given in Eq. (124) . To illustrate the 
essential features of the dynamics, in Fig. 3 we address the 
case of xq = 1 and show several snapshots of the atomic and 
photonic excitation probalities, i.e. respectively |(0| b x |*P(f))l 2 
and | <0| a x ^(f))! 2 , in the time interval te [0, 4tt/£2] (afterwards 
the same behavior is cyclically re-exhibited). The exciton ini- 
tially present at cavity 1 (so that at t = only the atomic bound 
mode is excited) progressively spreads over nearby odd cav- 
ities in the form of both field and atomic excitations. Such 
a stage, during which the field bound mode absorbs energy 
from the atomic one, finishes at time t-n/Q. when the field 
bound mode attains its maximum amplitude. Next, the field 
mode returns energy while the exciton probability spreading 
continues. At t — 2n/Q., the field is again fully unexcited 
but, remarkably, the exciton is no more localized on the first 
atom. Rather, for each atom in an odd cavity (except cav- 
ity 1) the excitation probability reaches its maximum value. 
Finally, between t-ln/Q. and t -An 1^1 while the field under- 
goes a further excitation-deexcitation cycle analogous to the 
previous one the excitonic distribution progressively localize 
around the first cavity until at t-An/Q. the initial state is fully 
retrieved. All of such features can be easily and accurately 
predicted once the moduli of Eqs. ( f30b and OTT l are squared, 
which when A = yields 



(2NJ\ 2 x+x2 . 2 £l 

PaM) = 



Sx,xo+Nt 1 I cos— t-l 



(32) 
(33) 



where pf x (p a ,x) stands for the photonic (atomic) probability 
excitation. 

It is worth mentioning that one can give a pictorial descrip- 
tion of the above dynamics in terms of breathing-mode be- 
haviors. While the field bound mode exhibits pure "transverse 
breathing" (with respect to the array axis), the excitonic mode 
in addition to this also shows "longitudinal breathing" since 
the atomic excitation somehow cyclically propagates along 
the array axis and localizes again on the starting atom. 

The setting r\ = -0.25 in Fig. 3 was chosen in order to bet- 
ter highlight these phenomena. Higher values of 77 reduce the 
number of involved cavity sites [we recall that the bound mode 
characteristic length is given by ^1=1/ log[( 1 +77)/ ( 1 - 77)], see 
Eq. ( fT3] l and Fig. 2]. For lower values of 77 more cavities get 
involved in the dynamics but the overlap between the exci- 
tonic bound mode and the initial state shrinks. As the compo- 
nent of the initial state that is orthogonal to the bound mode 



remains frozen in the light of Eq. (l24l . such overlap clearly 
affects the maximum amount of energy that can be exchanged 
between atoms and photons. This can be seen in Fig. 4, where 
for different 77' s we plot the overall atomic (photonic) excita- 
tion probability J^ x p a ,x (T*xPf,x) against time. Remarkably, 
in full analogy with a standard JC model IToTl the atom-field 
exchange of energy occurs so that the initially excited atom 
fully retrieves the energy released to the array regardless of 
the array size. Notice that this takes place at a rate given by 
the very same Rabi frequency associated with a single isolated 
cavity. As anticipated, the higher 77 the larger is the exchanged 
amount of energy. In particular, in the limit 77 — » — 1 the field 
is able to absorb the entire energy initially stored in atom 1. 
In such a limit, indeed, the coupling between cavities 1 and 
2 is fully broken. Hence, a standard JC dynamics takes place 
within cavity 1 . This suggests that one can regard the setting 
presented in this work as a sort of degenerate single-cavity JC 
model, which can describe the intermediate situation between 
a standard single-cavity JC and a uniform-hopping array. 

When the initial localized excitation is purely photonic, en- 
ergy exchange proceeds analogously. However, while the ex- 
citonic distribution evolves with time exactly as the photonic 
one in the above case (thus exhibiting pure transverse breath- 
ing) an analogous argument does not hold for the photonic dis- 
tribution. This exhibits a behavior more complex than those 
in Fig. 3 and, remarkably, it spreads with time over the entire 
array. This can be seen with the help of Eq. (l24l : while the un- 
bound atomic normal modes all have the same frequencies, so 
do not the photonic ones. Hence, unlike the process in Fig. 3 
where the part of the initial wavefunction not overlapping the 
atomic local mode remains frozen, when the initial excitation 
is photonic the non-overlapping part freely propagates along 
the array in the form of photonic excitation. 



V. CONCLUSIONS 

In this work, we have considered an array of coupled cav- 
ities and atoms with interspersed hopping strengths. We 
have first presented analytical solutions for the normal eigen- 
modes of the free field Hamiltonian under open BCs to high- 
light the emergence of two continuous bands with a bound 
mode occurring at the center of their energy gap. In con- 
trast to uniform- and parabolic-couplings arrays ]8|] where in 
the strong-hopping regime the atomic dynamics is frozen, we 
have shown that depending on which cavity is initially excited 
a significant exchange of energy between atoms and photons 
can arise. The associated dynamics is basically the one oc- 
curring with a standard JC model, where the aforementioned 
photonic bound mode and its excitonic analogue play the roles 
of the cavity mode and two-level atom, respectively. Remark- 
ably, the Rabi frequency associated with such effective JC- 
type dynamics is the same as the one associated with a single 
isolated cavity. In real space, an excitation initially localized 
within one cavity periodically spreads over nearby cavity sites 
of the same parity, in the form of both photonic and atomic 
excitations, and then localizes back on the starting site so that 
the initial conditions are retrieved. Interestingly, there is an 
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FIG. 3: (Color online) Snapshots of the atomic (red dots) and photonic (blue dots) excitation probabilities against the cavity number at time 
instants within the range te [0, 4n/£l] and for the initial state where atom 1 is fully excited. We have set 77 = -0.25, N = 101 and (in units of 
J) k = 100, ojf = 1000, A = 0. Although the plots were obtained from numerical solutions of the exact Hamiltonian H, in practice they remain 
identical when Eqs. i3Q\ and J3U are employed. Only cavity numbers such that both excitation probabilities significantly differ from zero are 
shown. 




time time 



FIG. 4: Overall atomic (solid line) and photonic (dashed line) exci- 
tation probabilities against time for rj = -0.125, -0.25, -0.5, 1. The 
initial state is the one where atom 1 is fully excited. We have set 
N = 101 and (in units of J) k= 100, w f = 1000, A = 0. 



intermediate time instant at which while the field is fully un- 
excited and the initially localized exciton is spread over the 
characteristic range of the bound mode. 

In this work we have restricted to the case of odd N. Indeed, 
in our staggered tight-binding model with open BCs the even 
and odd cases cannot be treated on the same footing like with 
cyclic conditions lfl9tl . In the even case, while major features 
such as the presence of a band gap with a discrete level at 
its center still hold the discrete level, when present, becomes 
twofold ll2lll . We have thus focused on the odd case merely for 
the sake of argument in order to better highlight the physical 
effects that we have presented. A comprehensive treatment of 
both cases will be the subject of a future publication ll22tl . 

Localized (bound) normal modes often occur in solid-state 
physics l23ll typically in the vicinity of localized defects or 
impurities that break the translational invariance of the host 
lattice. In a similar vein, they also appear in various CCAs 



scenarios such as arrays with one |[l2tl or two lfl3|| impurity 
atoms, T-type arrays with a single impurity atom Ill4ll and in 
CCAs with one or two detuned cavities 111511 . Here, the bound 
mode responsible for the phenomena that we have presented 
arises in a somewhat different way since no impurities or de- 
fects are present. Rather, its emergence is essentially a pure 
boundary effect stemming solely from the finiteness of the ar- 
ray length (we recall that under cyclic BCs this mode is absent 
IU9I0 . As such, aside from the specific context here addressed 
the present work provides a paradigmatic example of bound- 
ary effects in a CCAs scenario. 

Concerning an experimental test of the phenomena pre- 
sented in this paper, arguments analogous to those recently 
discussed elsewhere for CCAs with controllable hopping 
strengths JH H3l hold here as well. It is important to point 
out that even though we have often assumed large-size arrays 
in our numerical examples all the discussed effects in fact do 
not depend on the number of cavities. This makes their obser- 
vation feasible even with a small-size array, a setting that is 
widely expected to become accessible in the imminent future. 

Finally, it is worth mentioning that in the one-excitation 
Hilbert space of the field the bound mode in Eq. ( fTTt in fact 
defines an invariant state of the free field Hamiltonian Hf. In 
particular, notice that it makes the hopping part of the free 
field Hamiltonian © effectively vanish (this is the reason why 
the Rabi frequency associated with the effective JC-type dy- 
namics is the same as the one associated with a single isolated 
cavity). Invariant subspaces are states able to inhibit the trans- 
port of excitations through a quantum network. As such, they 
play a major role, although indirect, in some models recently 
proposed to explain the observed highly efficient excitation 
transfer in light-harvesting complexes 12511 . Our findings pro- 
vide a novel mechanism and context (of course well different 
from the aforementioned ones) where invariant subspaces can 
also play a significant role in excitation transport. 
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where we have used (IA6I ) to replace the cosine sums. 

Eq. (IA21 > is thus proven (note that the sum remains un- 
changed if the upper bound (N+l)/2 is replaced with (N-l)/2 
since the term corresponding to m = (N+l)/2 clearly vanishes). 

To demonstrate (IA31 l. we first use Eq. ( IA5b in the case m ± 
in order to derive the identity 



Appendix A: Useful formulas 

In this Appendix, we prove the three identities (to be used 
in the following Appendixes) 

sin (fa: + & k ) = { s in[fe(jc-l)] - rsin(fa)} , (Al) 

N + 1 

y ^ sin(fa) sin(£ x) — — - — Skk' , 



(A2) 



j j jsin[£(x- 1)] sm(k'x) 



N+l 



+ sin[fc'(x- 1)] sin(fa)) = — - — coskSkt , (A3) 

where k = 2mnl(N + 1) with m = 1, (N— l)/2 (and simi- 
larly for k' with m! being the associated integer). Eq. ( lAlb 
is straightforwardly checked by expressing the sine function 
in terms of complex exponentials and then replacing e ±l ® k 
through Eq. ( fTol l as 



k ( 1 - ?]) e ikx (e- ik - r) - e - ikx (e ik - r) 



sin (kx + §k) = 



K(\-rf) 

£k 



2i 



sin[A:(x-l)] - Tsin(fa)} . (A4) 



As for the remaining identities, using the well-known sum for- 
mula for geometric series it turns out that 



.1=1 x=\ 



m = 

2(1 -cos k) 



(A5) 



By taking the real part of Eq. ( lA51 l the corresponding formula 
for the cosine sum is obtained as 



^ cos [(k+k)x-k] = Re 



\k±k')x 



= [-1 +(-l)" ,+m ']Re 

[-1 + (-l)"' + " / ] (cos fa-cos £') 
(cos fa-cos k') 2 + (sin k+ sin k') 2 
which holds for k + k' . For k-k' we obtain that 



,(A8) 



Yco S [(k±k)x-k]\ k= , = £i 1 cosW2x-l)] = (A9) 
^ ^^-cosfc 

x= 1 v 2 

In the upper case (corresponding to the case k + k') the sum 
vanishes since clearly so does the sum over both sin(2fa) and 
cos(2fa). Use of prosthaphaeresis formulas allows to decom- 
pose each product of sines in Eq. ( IA3b in terms of a cosine 
sum. Thereby, the first sum in the left-hand side of ( |A3t can 
be decomposed as 



V • nt ivi • %> cos[(k-k')x-k]-cos[(k+k')x- k] 

y sin[fc(x-l)]sin£ x— ,(A10) 

x=\ 2 

while the decomposition of the second sum in ( IA3I ) is obtained 
from ( lAlOb by exchanging k with k! . It is now clear that for 
any k + k! the left-hand side of dA3t vanishes given that an 
exchange of k with k' evidently transforms dA8b into its oppo- 
site. On the other hand, it is immediate to see that when k = k! 
Eqs. dA9l ) and flAlOt entail that the left-hand side of ( |A3t re- 
duces to (N+l)/2 cos k5k,k>- The proof of identity ( IA3b is thus 
complete. 



Appendix B: Orthonormality conditions 



2^ cos (fa) = -r +( -ir 



(A6) 



Upon application of prosthaphaeresis formulas each product 
of sines appearing in the sum ( IA2t can be decomposed as the 
sum of two cosines so as to yield that 

N±l N±l 

Y^i ^ cos[(k-k')x]-cos[(k+k')x] 
2^ sin(fa) sm(k x) = ^ ~ 

x=l x=l 

1(AM_[_ 1+( _1)«]} = AM m = m > 

m + m! 

N+ 1 



(A7) 



Here, we demonstrate that the set of discrete functions spec- 
ifying the normal operators in Eqs. ( fT3l ) and ( fT31 ) fulfill or- 
thonormality conditions. To this aim, we first set a suitable no- 
tation (to be used in the next appendixes as well). Firstly, we 
relabel each continuous-band normal mode associated with 
{k, +} as {k, ±1), i.e. we replace ± with the numerical index 
fi = ± 1 . The expansion of each normal annihilation operator is 
rewritten as 
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^ <Pi,x K = ^ <#\2je-1 «2x-1 + ^ ^,2.t a 2x , (Bl) 



A = 1 X— 1 JT= 1 

AT+1 AM 

A = 1 A - 1 A - 1 
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where by virtue of Eqs. ( fTST l and ([T5T l the discrete real func- 
tions tp's are defined as 



<P(,2x-l = Rt x 



<Pe,2x=0 



(B3) 



<Pkn,2x-i = /J-!Bsin(kx+&k), fk^ix =S sin (far) , (B4) 
where for compactness of notation we have set 



M = 



2 n 



rj-1 yr m -i ' 



N + 1 



(B5) 
(B6) 



We recall that t is defined according to Eq. dlOl) . Our goal is 
to prove that the set of AT discrete functions [<p^ x > <Pkfi,x) satisfy 
the orthonormality conditions 

N 

YjVCsW,* = 1 . CB7) 

.1-1 

2V 

^ Vf,.v 9V-t = , (B8) 
.1=1 

;V 

fkfi,x Veii'.x = S w ^ m ' ■ (B9) 



v=l 



Upon use of Eqs. (|B3t and dB5t along with the well-known 
sum formula for geometric series the following identity holds 



V T 2(^-1) = I Z. 

Zj T r 2 _ 1 



4'/ 



■5T 2 . (BIO) 



Hence, it turns out that 

N+l N±± 



v=l 



which proves Eq. dB7l ). 

Using now identity ( 1AU along with Eqs. dB3t and ( lB4t , the 
left-hand side of Eq. dB8b can be arranged as 



x=\ 1=1 

oc {r^'sinf^x- 1)] -T x sm(kx)} 

A/-1 A/+1 

oc JVsin(fct) - ^Vsin(fa) = . (B12) 



where we have used the fact that in the second sum on the 
last line either of the terms corresponding to x = 1 and x = 
(N+l)/2 vanishes [we recall that k = 2mn/(N+l), see Section 
Hill . Hence, Eq. ( IB8I 1 is proven. 



As for identity dB9K using 
Eq. dB9l ) can be written as 



the left-hand side of 



z 



y>kn,x<Pk'fi',x = ' 



fufu' ^ sin (kx+&k) sin (k'x+$k>) 
x=i 



£ 

+ ^ sin (kx) sin (k'x) 

x=l 



(B13) 



Upon use of Eqs. dT6b . (IA2t and dA3b , the first sum within 
square brackets on the right-hand side is given by 



2^ sm (kx+'&k) sm (k'x+§k') 



N+l 



(1+t -2-rcosfa 



where we have used the identity 



1+t -2rcosfc 



K 2 (l -7]) 2 



N+l 

Sw = — r- <***>, (B 14) 



(B15) 



which can be checked straightforwardly through definitions 
( TTUb and (fTTT i and application of half-angle formulae. Using 
Eq. dA2b along with Eq. ( IB14b . Eq. ( IB 13b takes the form 



<Pkti,x<Pk'n>,x=W ~ + S M S l¥ 

A=l 



(B16) 



which shows that identity dB9l ) holds. 

All of the three orthonormality identities dB7b - dB9t are thus 
proved. 



Appendix C: Normal modes in the special case ij = 

In this Appendix, we prove that the set of normal frequen- 
cies and corresponding annihilation operators appearing in the 
decomposition of Hj d"T2l) and defined in Eqs. ([T3])-(fT5]) re- 
duce to (0 and ([8]), respectively, in the special case rj — 0, as 
expected. 

Using Eqs. ( fTTT i and (fl4l . in the case that rj — the normal 
frequencies associated with and a^- become, respectively 



o>k+\ifO — w/+2kcos(£/2) 



(CI) 



t^kArfO — (Of—2iccos(k/2) — w/ + 2/ccos ^— — nj . (C2) 

According to identity dAU . proven in Appendix lAl for rj = 
the quantity sin(far + j?a) becomes proportional to the sum of 
two sines [r — > - 1 according to Eq. ( fTOt l . Upon application of 
the corresponding prosthaphaeresis formula we obtain that 



sm(kx+& k )\^o = 



2 sin (*£=*+*£) cos | 



sin 



2 cos | 



--ff|(2x-l) 



5 (2x-D 



(C3) 
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Also, sin(fct) can be arranged in either of the equivalent forms As for the bound mode ( fT3b . for 77 = this reduces to 



sin(kx) = sin 





k 






sin 


5 <2*)_ 


= sin 


IH (2x) 



(C4) 



In the light of Eqs. (fT5T l. dC3b and dC4b can be arranged as 



— X 

2 



at- 



V N+l t-l 

X- 1 



(C5) 

a* ■ (C6) 



Given that /< = 2jmi/{N + 1) with m = 1, (N-l)/2 (cf. Sec- 
tionHITb the normal frequencies (ICU and associated operators 
( |C5t coincide with the first (AT— 1)/2 operators and respective 
frequencies in Eqs. (0 and ®, respectively. 

As for the normal modes specified by Eqs. (IC2b and dC6b . 
one can replace the argument of the cosine and sine functions 
with its opposite. This is possible since while the cosine in 
Eq. ( lC2t remains unchanged the normal operator ( IC6b is mul- 
tiplied by an irrelevant global phase factor. The associated 
wave vectors thus become 



<k \ (n(N+l-m) 
12-T N+l 



nN n[(N+l)/2+l] 



N+l 



N+l 



(C7) 



Eq. ( IC7l i shows that the normal modes corresponding to 
Eqs. dC2b and ( IC6b coincide with the last (N-l)/2 modes 
specified by Eqs. © and ® . 



ael^o = lim 



iV+1 . 



77- 1 y r 

~2~ N 

> sin 

N+l ^-f 

X=l 



x=l 

n[(N+l)/2] 



N+l 



a x , (C8) 



thereby coinciding with © for m-{N +1)12 (up to an irrele- 
vant global phase factor). On the other hand, for m-(N+l)/2 
© trivially reduces to a>f. Hence, in the special case 77 = 
all of the normal modes specified by Eqs. <fT~3T> - Cfl~5T> correctly 
reduce to the normal modes with uniform hopping rates. 



Appendix D: Free-field normal-mode decomposition 

In this Appendix, we will explicitly check the validity of 
Eq. ( TT2l . i.e. the normal-mode decomposition of the free field 
Hamiltonian Hj. Due to the orthonormality conditions ( IB7b - 
(|B9t , the (real) matrix of the coefficients that define the normal 
operators [cf. Eqs. JB3fc -( IB4b l is orthogonal. Its inverse thus 
coincides with its transpose, which yields that 



(Dl) 



By using this equation to express operators {a x } in terms of 
normal operators, the free field Hamiltonian © becomes 



H 



f = <»f 



2 fix + Z Z I Z Wv^P*** + Z Z ^>*^-* I + h - c -) I 

Vx=l / kfi k'n' \x=l ) kn \x=l 

(N-l \ (N-l \ (N-l 



+2 



V A= 1 



t k[i k'fi' \x=l 



kft \ a- 1 



(ajfifc + h.c.) , (D2) 



where we have set 



facts, Eq. (ID2I) considerably simplifies as 



-*[l-(-l)^] 

<Pkn,x + <Pe,x Vfy,x+i 



fjjf x - fk t i,x+lfk'fi',x + fk t i,xlfk' t i',x+\, 



gkn,x 



(D3) 

(D4) 
(D5) 



The three sums over x between brackets on the first line of 
Eq. ( ID21 > coincide with the left-hand sides of the orthonormal- 
ity identities (lB7t-dB9]|. Moreover, due to Eq. (IB3T ) the coeffi- 
cient of al&t on the second line clearly vanishes. Using these 



H 



f = w / 



(N-l 



l & c+j] & l< + ZZ Z^fi* 



kfi k'fi' 



V-v=l 



(N-l 



kfi V X— 1 j 



a k^k'n' 

(D6) 



A comparison between Eqs. (ID6t and (fT2l with the help of 
Eqs. (fTTT i and ( TBl i shows that the proof of ( TTZb is now reduced 
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to demonstrating the identities 

N-l 

^Pxgk^x = , 
x=l 

N-l 



(D7) 
(D8) 



In addition, we need to prove that Sk is given by Eq. (fTTI) . 

As for the proof of ( 1D7| |. we use Eqs. ( IB3b and ( IB4l i to 
derive the following identities 



N- 1 

^ (± Vftpkn, x +\ ipe, x 

x=l 


N-l 
2 

= ± ^ <Pkji,2x<Pe,2x-\ ~ 
x=\ 


+ - , (D9) 

T 


N-l 

/ l (±l') X <Pkii,x<Pt,x+l 
.i=l 


N- 1 

2 

= ± / \ (Pfy,2(x-l)<Pt,2x- 
x=l 


= y, (D,10) 


where 






r = 


N-l 
2 

x=l 


(Dll) 


These identities alonj 


; with Eqs. (ID31> and ( 


D3Tl allow to 



straightforwardly arrange the left-hand side of ( ID7t in the 
form 



N-l /j x 

JJp*,?*^ = |- + 1 + -- ? 7jr= 



T 



y = , 



(D12) 

where we have used definition (TTOb . Eq. ( ID7b is thus demon- 
strated. 

To prove Eq. dP81 l. we consider the following identities 
holding for any x=l, (iV- 1)/2 

P2x-i<Pku,2x<Pk'ti',2x-i = -p'S 2 K sin(fex) (1 + 77) sin(/c'x+iVX (D13) 
P2x<Pkn,2x<Pk>»;2x+i = ~h'S 2 k sm{kx) (1 -77) sin[/c'(x+ !)+#*'], (D14) 



where we have used Eqs. (IB4I ) and dD3l ). Eqs. ( ID 131 ) and 
dD14b sum to 

P2x-Wkn,2x9k'n',2x-1 +P2x<Pkn,2x<Pk'ii',2»l = -fJ.'S 2 K Sm(kx)S , (D15) 

where 

S = (l+r})sm(k'x+# k ) + (l-7j)sm[k'(x+l)+&kl ■ (D16) 



Using now identities dAll i and dB15b along with definition 
( [Tol l, >S can be arranged in the form 

S= ^-^{[(l-J7)-T(l+^)]sin(fe'x) 



+(77+ 1) (sin[/c'(x- 1)] + sin[k'(x+ 1)])} 

£7 



: [( 1 _ T 2 _ 2r cos /c') sin(/t'jc)l = — sin(/cV) . (D 1 7) 

Sk' 1 J KEk' 



In deriving Eq. dD17| ) we have made use of prosthaphaeresis 
formulas to replace the sum of sines on the second line of 



( ID17b with the product 2cos k'sin(k'x). Using this result, upon 
sum of dD15t over x=\, (N— 1)/2 we thus obtain that 



Z 



jj.' ' Sk 5kk' 

P2x-lfkii,2xfk'ii',2x-l + p2xfkn,2xfk' n' ,2x+l = « , (D18) 



where we have used identity (IA21 > to carry out the sum of 
sin(/cx) sin(/c'x). 

Consider now the left-hand side of Eq. dD8b . where the sum 
can be split into the sums over even and odd terms. A com- 
parison of this with the left-hand side of Eq. (ID 18b should 
make clear that (|D8t can be obtained by adding dD18l ) to the 
quantity obtained from dD18b through exchange of {/c,yu} with 
{k' ,n'}. This argument leads to 



ZM _ )e k d k k> 
PxJkpi,x - 9 

x=l Z 



-p-Ekdkk'S, 



(D19) 



which proves that Eq. dD8b holds. 

Although we have proven both the identities dP7| ) and dP8l ), 
to finalize the explicit demonstration of Eq. (TTZb we need to 
prove that Sk is given by Eq. ( fTTI ) (note that in this Appendix 
up to this stage we have left Sk unspecified). To accomplish 
this task, we point out that Eq. ( fT6l i entails the constraint 



1 



K 2 (1-TJ) 2 



(P20) 



i.e. the squared modulus of the complex exponential on the 
left-hand side must be unitary. By bringing s 2 to the left-hand 
side and calculating the squared modulus on the right-hand 
side Eq. dP20b takes the form 



s \ = k 2 {[(77- l)cos/c- (77+ 1)] 2 + (t/-1) 2 sin 2 k] 



ok j . j k 
4/r cos" - + 77 sin - 
2 ' 2 



(P21) 



where in the last step we have made use of half-angle formu- 
lae. By taking the square root of dP21b Eq. (fTTI) is retrieved. 
The proof of decomposition (fT2l is therefore complete. 



Appendix E: Full Hamiltonian normal-mode decomposition 

Here, we prove that the full Hamiltonian (Q]) can be de- 
composed in terms of normal operators according to Eq. d20b . 
Having proven in Appendix [P] that (fT2l holds, it suffices to 
demonstrate that 



k M 



(El) 



x=l 



kfi 



x=l 



We first note that due to Eqs. ( fT8l and (fT9b the orthonormality 
conditions (IB7b - (|B91 > entail that Eq. ( IDlb holds for the atomic 
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operators as well, i.e. 



kji 



(E3) 



Using this equation to express each atomic site operator in 
Eq. (lElb in terms of normal operators along with the or- 
thornormality identities (|B7| i- (|B9l l, the proof of Eq. ( IE 1 b is 
in fact identical to the one carried out in the case of field op- 
erators [cf. first line of Eq. ( ID2l i and the quantity between the 
first brackets on the right-hand side of Eq. flD6H . Likewise, 
use of Eqs. (|B7]i-(|B9]l and (|E3j yields that 



By summing Eq. (IE4t to its adjoint Eq. (IE2b is retrieved. The 
proof of normal-mode decomposition ( f2Qb is thus complete. 
Note that the uniformity of the atomic frequencies and atom- 
photon coupling strengths throughout the cavity array is cru- 
cial for Eqs. (IEU and ( IE2b . and hence d20l i. to hold. 



x=\ 
( N 



x=\ x=l \ kii )\ k'u' 



Zj&x & \Pe + ZaZu J_ I <Pk f i,x<Pk> fi ',x 

(N \ 
.1=1 / 



(E4) 
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